Trap anharmonicity and sloshing mode of a Fermi gas 
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For a gas trapped in a harmonic potential, the sloshing (or Kohn) mode is undamped and its 
frequency coincides with the trap frequency, independently of the statistics, interaction and tem- 
perature of the gas. However, experimental trap potentials have usually Gaussian shape and anhar- 
monicity effects appear as the temperature and, in the case of Fermions, the filling of the trap are 
increased. We study the sloshing mode of a degenerate Fermi gas in an anharmonic trap within the 
Boltzmann equation, including in-medium effects in both the transport and collision terms. The 
calculated frequency shifts and damping rates of the sloshing mode due to the trap anharmonicity 
are in satisfactory agreement with the available experimental data. We also discuss higher-order 
dipole, octupole, and bending modes and show that the damping of the sloshing mode is caused by 
its coupling to these modes. 

PACS numbers: 67.85.Lm 



I. INTRODUCTION 

In a couple of experiments, the measurement of the 
frequencies of collective modes in trapped Fermi gases re- 
vealed a lot of interesting information on the equation of 
state, the validity of superfluid hydrodynamics, and the 
superfluid- normal phase transition [lH5j- More recently, 
the transition from the hydrodynamic to the collision- 
less regime in the normal phase was also studied d, 0] ■ 
However, since the frequencies of the collective modes de- 
pend on the trap frequencies, a precise knowledge of the 
latter is required for a meaningful interpretation of the 
collective-mode data. 

A possibility to determine the trap frequency with high 
precision is the measurement of the frequency of the 
sloshing mode which is an oscillation of the center 
of mass of the trapped atom cloud. In a harmonic trap, 
this oscillation (also called Kohn mode) is undamped and 
its frequency coincides exactly with the trap frequency in 
the corresponding direction, independently of the num- 
ber of atoms N, or their interaction, their temperature T 
etc. This is a consequence of the Kohn theorem and 
follows from the fact that the center-of-mass oscillation 
decouples completely from the internal dynamics of the 
gas if the interaction is translationally invariant and the 
external potential is harmonic |10| . 

However, in practice the trap potential is never exactly 
harmonic. In optical dipole traps the potential is 
typically Gaussian (corresponding to the intensity profile 
of the laser beam). Because of this anharmonicity, the 
frequency of the sloshing mode is shifted, and the shift 
depends on the system parameters such as N, T, the scat- 
tering length a characterizing the interaction strength, 
etc. Furthermore, the sloshing mode is no longer un- 



damped. 

The anharmonicity of the trap complicates consider- 
ably the analysis of collective-mode experiments. In Ref. 
Q the measured frequencies were corrected for the an- 
harmonicity effects by giving them in units of the mea- 
sured frequency of the sloshing mode. However, it is 
clear that the damping rates of the modes cannot be cor- 
rected in this way. In Ref. [l2j], the damping rate of the 
sloshing mode was used to estimate the increase of the 
damping rate of other modes due to the anharmonicity. 
Both corrections are ad-hoc prescriptions without rigor- 
ous justification. It is therefore strongly desirable to get 
a better understanding of the anharmonicity effects on 
the sloshing mode. 

The aim of the present paper is to describe theoret- 
ically the frequency shift and the damping rate of the 
sloshing mode in an anharmonic trap. We will compare 
our results with the experimental data available from the 
Innsbruck group @, 03 and with the numerical results 
by Wu and Zhang We will also discuss in detail 

the damping mechanism of the sloshing mode. In an 
anharmonic potential, the center-of-mass motion is no 
longer decoupled from the internal degrees of freedom of 
the cloud. We will see that the damping of the trans- 
verse sloshing mode is a consequence of its coupling to 
other damped collective modes, in particular to the radial 
dipole mode and the bending mode. 

The framework of our study is the Boltzmann equa- 
tion, including mean field [15| and in-medium cross- 
section 

El. Especially the mean field is expected 
to be important in the present context, because it can 
have a sizable effect on the density profile, i.e., on how 
far the cloud extends into the anharmonic region of the 
trap potential. The Boltzmann equation is solved ap- 
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proximately with the help of the phase-space moments 
method. This method, when extended beyond the lowest 
order, has proven to be in very good quantitative agree- 
ment with the results of a full numerical simulation |17| . 
It has also been quite successful for the description of the 
experimental results for the fre quen cy and damping rate 
of the radial quadrupole mode [18| . In the present case 
of the sloshing mode, we include phase-space moments of 
first and third order. 

The paper is organized as follows. The general formal- 
ism is briefly presented in Sec. [TTJ Then we specify our 
model for the experimental trap potential in Sec. IIII1 In 
Sec. IIV1 we give a formula for the frequency shift of the 
sloshing mode within the first-order moments method. 
Then we extend the ansatz to third order in Sec. [V] The 
physical contents of the extended ansatz and the numeri- 
cal results are discussed. In Sec. lVII the third-order ansatz 
is used to describe also the radial dipole, radial octupole, 
and bending modes. Finally, in Sec. IVII1 we will con- 
clude. 

Throughout the paper, we use units with h = ks = 1- 



II. SUMMARY OF THE FORMALISM 

In this section, we give a short summary of the for- 
malisms of Refs. [ill [l8j]. More details can be found 
there. 



Linearized Boltzmann equation with in-medium 

effects 



We consider a balanced two-component (Nf — N± = 
N/2) Fermi gas of atoms with mass m and interspecies 
attractive interaction (scattering length a < 0), trapped 
in a potential Vt{y). The framework we use to describe 
the collective dynamics of the system in the normal-fluid 
phase is the Boltzmann equation. We assume that the 
two components move in phase, so that only one distri- 
bution function / = /-j- = /^ is needed. It is normalized 
to J d 3 rd 3 p/(2ir) 3 f — N/2, and expectation values of 
one-body operators are given by 



(9>(*) 



2 

iV 



d 3 rd 3 p 

-^jr/(r ) p,t)«(r,p). 



(1) 



As in Ref. [la ] ■ we include in- medium effects in both the 
transport and the collision parts of the Boltzmann equa- 
tion: a mean field like potential U and the in-medium 
modified cross-section da/dVL, respectively. 

Within the Thomas-Fermi or local-density approxima- 
tion (LDA), the equilibrium distribution function reads 



/«/(*■» P) 



1 



(2) 



where U eq is the mean field in equilibrium, /iq is the 
chemical potential, and /3 = 1/T is the inverse temper- 



ature. As in Ref. [l5|, we obtain U eq from the single- 
particle self-energy in ladder approximation, evaluated 
at the Fermi level. In the weak coupling limit, this re- 
duces to the Hartree term t/Hartrcc = 47rap/m, where 
p = J d 3 p/{2-jr) 3 f is the density per spin state, and it 
remains finite for all interaction strengths up to the uni- 
tary limit, a — > — oo. In equilibrium, the main effect of 
the mean field is to enhance the de nsity in the center of 
the trap, as shown in Fig. 3 of Ref. [l5| . 

For the study of collective oscillations, it is sufficient 
to consider small deviations from equilibrium and to lin- 
earize the Boltzmann equation with respect to Sf — 
f ~ feq- If we write the variation of the distribution 
function in the form Oil 



Sf(r,p,t) = / eg / eg $(r,p,£) . 



(3) 



with / eg = 1 — / e g, the linearized Boltzmann equation 
can be written as 



/eg/eg ^+{^^-+V T + U eq } 



V r (SV + SU) 



-/[$], (4) 



where {F, G} = V r F ■ V P G - V P F ■ V r G denotes the 
Poisson bracket, SV is the perturbation of the trap po- 
tential that is used to excite the collective mode, SU is 
the variation of the mean field due to the variation of the 
density, and J[$] is the linearized collision term. Since 
we want to calculate the so-called response function, we 
take the perturbation to be a pulse, 



5V(r,t) = V(r)S(t). 



(5) 



As in Ref. [l5j , we approximate the variation of the mean 
field by 



SU(r,t) 



dp 



Mr,*). 



(6) 



Pe,(r),T 



The linearized collision integral reads 



d 3 pi 
(2tt) 3 



do- |p- Pi| 



feqfeqlfeqfeql 



dfl m 

X + (7) 



where an obvious notation for the different f eq and $ 
at the momenta before (p, p x ) and after the collision 
(p'j Pi) nas been used. Note that, especially near the 
critical temperature, the in- medium cross-section da/dfl 
can differ strongly from the free one @, [H, 13 • 



B. Moments method 

As in Ref. [l8[, we are looking for a semi- analytical so- 
lution of the Boltzmann equation Q by using the method 
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of phase-space moments. In this section, we will general- 
ize the formalism of that paper to the case of an arbitrary 
trap potential Vt and with mean field U. 

The basic idea is to approximate the function $ by 
a polynomial in the components of r and p with time- 
dependent coefficients Cj, 



$(r»P,*) = ^c i (t)0 j (r,p) 



(8) 



where the <pj are suitable basis functions, e.g., monomials 
in the components of r and p. 

Multiplying the linearized Boltzmann equation (|4]) by 
(j>i and integrating over phase space, one obtains, after a 
Fourier transform with respect to t, a set of n coupled 
linear algebraic equations for the n coefficients Cj(uS). In 
matrix form, they read 



(9) 



where 



- —iujMij - 


i a trans i aSU , a coll 
t- Jiij + /iy -t- Ajj , 


(10) 


Mij = J 


d 3 rd 3 p f - x x 
(27r) 3 eq * ' 


(11) 


a trans 


N 

~lj»({<l>i,<t>j})eq , 


(12) 


A su_N 


(Vp4>i ■ V r SUj) eq , 


(13) 


A COll 

A ij ~ 


f d 3 rd 3 p 

1 (2*)» * W ' 


(14) 



and 



V r F(r)) { 



(15) 



The variation of the mean field in Eq. (1131) is defined as 

d 3 p 



dp eq J (2?r) 



rfeqfe 



(16) 



Equations (TT21 and (IT5|) correspond to Eqs. (14) and (16) 
in Ref. [18] which have been simplified by integration by 
parts. Comparing Eqs. (|9|- (fT5l) with the analogous ones 
of Ref. [Hj], Eqs. (11)-(16), one sees that the mean field 
gives rise to a new term, A 5U . However, one should keep 
in mind that implicitly all terms depend on the mean 
field since it modifies the equilibrium distribution f eq . 
One can see that the matrices M and A co11 are symmet- 
ric, A trans is antisymmetric, whereas A su has no defined 
symmetry. In practice, the calculation of the matrices 
is straight-forward, but tedious, and we made use of the 
Mathematica software to express the numerous matrix 
elements in terms of a smaller number of integrals over 
equilibrium quantities. 



C. Eigenmodes and response function 

Without an external perturbation V, i.e., for a, = 0, 
Eq. j9|) has a solution with non- vanishing coefficients Cj 
only if det A(lo) = 0. The frequencies lo for which this 
happens are obviously given by the eigenvalues of the 
matrix -iM~ l (A trans + A su + A coU ). If they are well 
separated, it is possible to interpret them as the frequen- 
cies of the eigenmodes of the system. In general, they are 
complex, and their imaginary part describes the damping 
rate of the corresponding mode M,ll5l ] . 

However, as discussed in Ref. fl8j |. when the moments 
method is extended to higher order, there can be many 
eigenvalues belonging to a single collective mode. In this 
case, the scattering of the eigenvalues, which goes over 
into a continuous spectrum in the limit of an infinite num- 
ber of moments [201 ], corresponds to a new contribution 
to the damping (Landau damping) in addition to the 
imaginary parts coming from the collision term. In order 
to obtain the mode frequency and damping rate in this 
case, it is useful to look at the response function which 
contains the contributions of all eigenvalues. 

We denote by (q)(u>) the Fourier transform of the ex- 
pectation value (q)(t) of some operator q after the per- 
turbation. The so-called response function is equal to 
S(q)(uj)/a = ((q)(uj) — (q) eq )/a in the special case that 
the excitation operator is V = aq. In all excitations con- 
sidered in this paper, (q) eq will be zero. The strength 
function is proportional to the imaginary part of the re- 
sponse function. 

In order to calculate the response function, we need 
the coefficients Cj(ui) of the ansatz ([8]). By diagonalizing 
the matrix 



\5U 



A co11 ) = PDP- 1 , (17) 



M~ 1 (A trans + A 
with D = diag(ri . . . , T n +iu> n ), we can write them 



as 



. P jk (p- 1 M- 1 a) k 



fc=i 



(18) 



Then, the response function can be easily obtained as 
(q)(u>) = b T c(cu), where 



N 



d 3 rd 3 p 
(2tt) 3 



feqf- 



qj eqH*rt - 



(19) 



The frequency and damping of the mode can be ex- 
tracted from the response (q)(u>), e.g., by fitting the peak 
in the strength function — lm(q)(uj) corresponding to the 
mode under consideration with a Lorentzian. 

Using the solution for the coefficients Cj (lo) at the peak, 
we can also obtain the velocity field of the corresponding 
collective mode: 



" 1 f d 3 p 

v(r >W )=5>(o,)^- J J^fe^ 



(20) 
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III. REALISTIC TRAP POTENTIAL 

In this section we recall the shape of a typical optical 
dipole trap. We concentrate on the configuration used 
by the Innsbruck group, consisting of a focused-beam 
trap [TTj] with additional magnetic confinement in the 
axial direction. Following Ref. [2lJ, we parametrize the 
experimental trap potential as 



V exp (r) = Vo 



1 



■ exp 




ho.x 



(21) 



where r± = \J x 2 + y 2 , Vo is the trap depth, wq is the 
minimal waist of the laser beam, zq — ttwq/X is the 
Rayleigh length, A the laser wavelength and uj z ,mag de- 
fines the magnetic trapping in the z direction. 

In order to compare our results with the data from 
Refs. 0, we focus on the setup of that experiment. 
The trap was of the type (|21j) . with strong anharmonicity 
in the x-y plane due to the Gaussian shape of the laser 
beams. The confinement along the z direction was prac- 
tically harmonic. In our calculations, we will expand Eq. 
(1211) up to second order in z and up to sixth order in r± . 
The result can be written as 



V 6th (r) 



mu 2 r 2 f m.j" r 



lr 2 



W 



9 4 4 

24y 2 



FIG. 1: Experimental trap potential (solid line), and its har- 
monic (dashes) and sixth-order approximation (dash-dotted 
line), as functions of x for y = z = 0. The Fermi energy for 
the parameters of Ref. 0] is indicated by the thin dashes. The 
potentials are in units of lo x and x is in units lho,x = l/y/rnuj x . 



As pointed out in Ref. pfl , the anharmonicity effects de- 
pend mainly upon the ratio Ef/Vq. Defining the Fermi 
energy as usual by Ep — (SiVor^Wz) 1 / 3 , one obtains 
Ef/Vo ~ 0.1, i.e., at low temperature the anharmonicity 
effects are relatively weak. For illustration, Ep is indi- 
cated in Fig. [T] by the thin dashed line. We see that, 
for the present choice of parameters, the atoms start to 
feel the anharmonicity of the potential when their energy 
exceeds Ef, while the sixth-order approximation to the 
potential stays very precise up to about five times the 
Fermi energy. 



(22) 



Wo/{mwl) and oj 2 = + 



z,opt 



where 

w 1,mag + ^■Vol{ mz o) are the radial and axial trap fre- 
quencies, respectively. We use Eq. (|22[) instead of Eq. 
(f2~Tj) because it simplifies the calculations in the sense 
that it leads to a density that tends to zero for r± — > oo. 
The necessity to go beyond fourth order is also clear from 
Eq. (|22|) : otherwise the potential would be unbound from 
below. For illustration, Fig. [1] shows the real trap poten- 
tial (|2ip . its harmonic approximation, and the one we 
will use, Eq. ([22j) . 

In the present paper, we concentrate on axially sym- 
metric traps. Nevertheless, we write u> x and oj y instead 
of lo±_ for the trap frequencies in x and y direction if the 
formulas can be generalized to the triaxial case. 

In our numerical calculations, we use as an example 
the parameters of the Innsbruck experiment Q. As 
mentioned in Ref. [7J, the sloshing mode was studied 
with the same parameters as the compression mode, i.e., 
cj ± /(2tt) = 1100 Hz, w z /(2tt) = 26 Hz, and Vo = 19 fM. 
The trap was loaded with N — 6 x 10 5 atoms of 6 Li in 
the unitary limit, l/(kpa) — 0, and the temperature was 
varied between ~ and \ flTp. Since we cannot describe 
the superfluid phase, we limit ourselves to temperatures 
above 0.3Tp [15| . Moreover, we approximate the unitary 
limit numerically by setting l/(kpa) = —0.01. 



IV. SLOSHING MODE AT FIRST ORDER 

In order to describe the sloshing motion, say, along 
the x direction, the ansatz (jSJ has to contain at least two 
basis functions: 



$j„t(r,p,*) = d(t)x + c 2 (t)p x 



(23) 



The second term describes the collective velocity of the 
cloud in x direction, while the first one corresponds ap- 
proximately to a displacement of the center of mass. The 
latter statement becomes exact in the special case of a 
purely harmonic trap without mean field (£7 = 0). In 
this case, the ansatz (1231) is closed with respect to the 
operator p/m • V r — V r Vr • V p on the left-hand side of 
the Boltzmann equation and since the collision term 
does not contribute (I[x] = I\p x ] = because of momen- 
tum conservation in a collision) , the ansatz (|23|) becomes 
exact. However, if the trap is not purely harmonic, this 
is no longer true, since the gradient of the trap potential 
generates new terms. And even if the trap is harmonic, it 
is not true for interacting atoms, which feel a mean field 
U ^ in addition to the trap potential Vr- In the next 
section, we will therefore extend the ansatz (|23[) . but let 
us first look what happens at first order. 

We now apply the formalism described in Sec. [TT] to 
&ist- The calculation of the matrix A is very easy in this 
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case. After some algebra, we obtain a quadratic equa- 
tion for the frequency of the sloshing mode by imposing 
dot A = 0. Its solution reads 



10 



1-C 



losh 



1-C-C 



The parameters 



and 

C ■ 



2 A 

' N 



c- 2 l 

N 



d 3 rd 3 p 



d 3 rd 3 p 
(2tt) 3 



feqft 



dU, 



eq 



qj eq-* 



dx 



(24) 



(25) 



2„,2 



(27T) 



feqfe 



d /miL'-.c 
dx \ 



Vt(t)) (26) 



characterize, respectively, the strength of the mean field 
and of the anharmonicity effects. Equation ([24| has two 
important features: (a) If Vr is harmonic (C = 0), then 
u slosh = w x , independently of the interaction, in accor- 
dance with the Kohn theorem Q . This point was already 
discussed in Ref. [HI and shows the consistency of our 
approach, in particular of Eq. ([6]). (b) If Vr is not purely 
harmonic (C ^ 0), then u> s iosh depends on both the an- 
harmonicity of the trap and on the interaction. 

Another property of Eq. (IM|) is that the sloshing mode 
is undamped (i.e., ui s i os h is real) even in the case of an an- 
harmonic potential. This is, however, only a consequence 
of the first-order approximation. 

A similar calculation was performed in the Appendix 
B of Ref. Q- In that reference, only the leading anhar- 
monicity correction was kept. If we combine Eqs. (|22[) . 
(|2~4")) . and (l2"rJ|) and keep only the leading order in 1/Vb, 
we find 



2 

w slosh 



muj]_(i{x 2 ) eq + (y 2 )eq) 

2V (1 - C) 



(27) 



This result differs from Eq. (B2) of Ref. [7] in two re- 
spects. First, there is an additional factor of three in front 
of the (x 2 ) term (misprint in Ref. [?})■ Second, the an- 
harmonicity correction is enhanced by a factor 1/(1 — C). 
This factor is missing in Ref. Q because there the mean 
field was not considered, although the expectation values 
(x 2 ) and (y 2 ) were calculated with the density profiles of 
an interacting gas. 



V. SLOSHING MODE AT THIRD ORDER 

A. Extended ansatz 

As in Ref. [lH, we will extend the ansatz (f2"3")l by in- 
cluding higher-order terms. At the next higher order, the 
ansatz contains 18 terms 

18 

$ Srd (r,p,t) = ^c i (^ i (r,p), (28) 
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FIG. 2: Strength function — Im(a;)(w) for the sloshing mode 
(excitation V = ax) in the anharmonic trap for various tem- 
peratures. The system parameters are those of Ref. 0. 



where: 



x , 02 = Px , 



x 2 p x , 



>5 = xp x , 6 



b 7 = xy 2 , 8 = y 2 p x , (j>g = xypy , <f> 10 



■ Pxpl , 013 = xz 2 
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pl, 

: VPxPy 
= z 2 p x , 



Pn = xPy , 9>i2 

f>15 = XZp z , 016 = Zp x p z , 017 = xp\ , 018 = PxPt- (29) 



The first two terms of $3 rt j coincide with $>xstj & n( i the 
subsequent ones are all possible terms of third order in 
the components of r and p satisfying the symmetries of 
the mode: odd under (x,p x ) — > (—x, —p x ) and even un- 
der (y,p y ) -> (~y,-p y ) and (z,p z ) -> (-z,-p z ). 

Note that some of the third-order terms, e.g., 04 = 
x 2 p x , generate a velocity field that is quadratic in the 
coordinates, whereas the velocity field corresponding to 
the first-order ansatz is spatially constant. Other terms, 
e.g., 05 = xp 2 or 06 = p 3 , describe momentum-sphere 
distortions. They give rise to a non- vanishing collision in- 
tegral so that they generate a damping. The momentum- 
independent third-order terms such as 03 = x 3 describe 
deformations of the cloud in coordinate space that will 
be discussed in more detail in the next subsection. 

As explained in Sec. Ill CI the frequency and damping 
rate should now be extracted from the response function. 
The suitable excitation operator is 



Vsiosh (r) = ax . 



(30) 



because it gives all atoms a constant kick in x direction, 
Px — >■ Px — a, at the moment of the excitation t = 0. 
The response function is defined as the expectation value 
(x) {uj)/a. 

We checked numerically that, in a purely harmonic po- 
tential (Vb — > oo), the strength function has only a single 
sharp peak (no damping) at u> = io x within the precision 
of our calculation (better than 10 -4 ), independently of 
the temperature and of the interaction strength. 

In Fig. [3 we show the strength functions for the slosh- 
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FIG. 3: (a) Temperature dependence of the frequency of the 
transverse sloshing mode in units of ui x . Solid line: third- 
order result; dashed line: first-order result; dash-dotted line: 
first-order result obtained without mean field. The experi- 
mental data points are taken from Ref. [7|. (b) Temperature 
dependence of the damping rate in units of ui x . Solid line: full 
third-order result; dashed line: third-order result obtained 
with the free instead of the in-medium cross section. The 
filled and empty data points are measured damping rates of 
the sloshing modes in x and y direction, respectively, from 
Ref. O. 



ing mode in the anharmonic trap for various tempera- 
tures. In all cases the strength is concentrated in a sin- 
gle peak and the width of the peak corresponds to the 
damping rate. The position of the peak is always be- 
low uj = lj x , and with increasing temperature, this shift 
gets stronger. This is easily understood: With increasing 
temperature, more and more atoms reach the peripheral 
region where the anharmonic potential is flatter than the 
harmonic one (cf. Fig. [1]) . 

Since in all the cases shown in Fig. [2] the strength is 
concentrated in a single peak, this allows us to extract 
the frequency and damping rate by fitting it with a sin- 
gle Lorentzian. (This is not clear a priori: For instance, 
if we do the calculation without mean field, the strength 
function for low temperature has two separate peaks, and 
it is not evident how one has to define the average fre- 
quency and total width.) The results for the frequency 
and damping rate as functions of temperature are shown 
in Fig. [3Ja) and (b) as the solid lines. For comparison, 
the dashed line in Fig. [21a) shows the first-order result 
for the frequency. The dash-dotted line corresponds to 



the first-order result obtained without mean field. As ex- 
pected, the mean field reduces the shift of the frequency, 
but its effect is moderate and most pronounced at lower 
temperature. 

The experimental results from Fig. 2 of Ref. 0] are also 
shown in Fig.^a). In that paper, the measured empirical 
temperatures T were already converted into real temper- 
atures T/Tp as explained in [H, [23[, so that the results 
can immediately be compared with our calculation. The 
trend is correctly reproduced, and the third-order result 
is clearly in better agreement with the data than the 
first-order one, but the theoretical frequency shift is still 
stronger than the experimental shifts reported in Ref. 0] ■ 
However, one should keep in mind that while Lo s iosh was 
precisely determined, the ratio uj s i os h/u x depends also on 
uj Xi which was not very well known in that experiment. 
Actually, cu x was deduced from u) s iosh under the assump- 
tion that the anharmonicity effects become negligible at 
T = 13]. If we suppose that the true trap frequency oj x 
was just ~ 1.5% higher, the data points for u s i os h/oj x are 
slightly shifted downwards and the agreement between 
theory and data becomes much better. 

The damping rate, shown in Fig. [3Jb) , increases 
strongly with the temperature: at T/Tp = 0.4 the slosh- 
ing mode survives for several hundreds of oscillation peri- 
ods, whereas at T/Tp = 1.2 its amplitude decreases by a 
factor of 1/e after ~ 6 oscillations. The in-medium modi- 
fication of the cross-section has only a small effect on the 
damping of the sloshing mode. The experimental data 
from Ref. [13| are in quite good agreement with our the- 
oretical result, especially at the high temperature. The 
strongly scattered damping rates at low temperature are 
probably an artefact due to a beat caused by the small 
residual ellipticity of the trap potential in the experiment. 

We also compared our results with those shown in Fig. 
5(a) of the recent paper [l4| by Wu and Zhang, where 
they solve the Boltzmann equation numerically within 
the relaxation-time approximation with a local relaxation 
time r(r) taken from Ref. [24(. As we discussed in our 
previous work [l?], Oil , the inclusion of higher-order mo- 
ments accounts in an approximate way for the spatial 
dependence of the relaxation time. For the sake of com- 
parison, we use the same trap parameters as Wu and 
Zhang 1 , namely, uj±/(2tt) = 1800 Hz, oj z /{2n) = 32 Hz, 
and Vq = 50 fjK, and neglect, as it is done in Ref. [l4j . 
the mean field and in-medium effects on the cross section. 
Because of the larger trap depth, the frequency shifts are 
somewhat weaker than those shown in Fig. [2£a). Since 
in the calculation by Wu and Zhang the frequency de- 



The comparison made in Fig. 5(a) of Ref. [3 with the data 
of Fig. 2 of Ref. Q is somewhat misleading because the trap 
parameters are different. In the calculation of Ref. [l4| . the pa- 
rameters corresponding to the quadrupole-mode measurement of 
Ref. Q were used, whereas in Ref. Q it is said that the sloshing 
mode was studied with the parameters of the compression-mode 
measurement. 



7 



(a) — 




FIG. 4: Schematic representation of the oscillation of the den- 
sity profile in x direction (p x being the density integrated over 
y and z) for (a) the Kohn mode, in which the cloud oscillates 
as a whole, and (b) the radial compressional dipole mode, 
where the cloud shape oscillates, while the center of mass of 
the cloud stays at rest. All density profiles have the same 
normalization (number of atoms) and the same width. 



(a) 



(b) 



FIG. 5: Schematic representation of (a) the Kohn mode, in 
which the cloud oscillates as a whole, and (b) the bending 
mode, where the cloud shape oscillates, while the center of 
mass of the cloud stays at rest. 



pends on the oscillation amplitude, we have to compare 
with the results obtained for the smallest amplitude [red 
squares in Fig. 5(a) of Ref. [13]] ■ If we extrapolate these 
results, shown only up to T/Tp = 0.4 in Ref. [1J|, to 
higher temperatures, they are in quite good agreement 
with our results above T/Tp ~ 0.5. At lower temper- 
atures, however, our frequency shift is slightly stronger 
than that by Wu and Zhang, probably because they use 
Boltzmann instead of Fermi distributions so that their 
density is more concentrated in the trap center at low 
temperatures. 



B. Coupling between sloshing and other modes 

In a harmonic trap, a characteristic feature of the Kohn 
mode @ is that, as long as the interaction is translatio- 
nally invariant, the cloud moves as a whole without any 
change in size or shape |io| . as illustrated in Figs. QJa) 
and[5ja). In an anharmonic trap, this is no longer the 
case and the cloud will be slightly distorted during the 
oscillation. In other words, because of the anharmonicity, 
the pure center-of-mass motion gets coupled to other col- 
lective modes. These are damped as usual by collisions, 
which results in a damping of the sloshing mode. 

What kind of distortions of the cloud shape can one 
expect? Let us first note that, to first order in the per- 
turbation, the cloud width cannot be changed (neither 
in the transverse nor in the axial direction), because the 



excitation operator is odd with respect toi<4 —x. 

However, the density profile can become skewed in the 
direction of the oscillation (i.e., in the x direction in our 
case) . For illustration, Fig. 0Jb) shows how the density 
profile changes during a radial compressional dipole oscil- 
lation. To derive the corresponding operator, let us start 
from the usual definition of the skewness, which is pro- 
portional to ((x— (x)) 3 }. Keeping only terms linear in the 
variations 5(x) and S(x 2 ), and using (x) eq = 0, we find 
that for small oscillations, this is equal to (x 3 — 3(x 2 ) eq x). 
If the trap is axially symmetric, collective modes of dif- 
ferent multipolarity in the radial direction do not mix. It 
is therefore useful to decompose the operator into dipolar 



and octupolar parts: a; 3 — 3(x 2 



(:!(/•' 



q oct with 



4 y oc 



qdip = {r± - 2{r 2 L ) eq )r 1 _ cosip = (x 2 + y 2 - 2(r 2 L ) eq )a 



q oc t = cos 3(p = (x 2 - 3y 2 )x . 



(31) 
(32) 



Only qdip can get a non- vanishing expectation value dur- 
ing the sloshing oscillation. 

Another possible distortion involves the axial (z) di- 
rection. Let us consider again the sloshing mode in x 
direction. Then one could imagine that the cloud near 
z = 0, owing to its larger radial size, is more sensitive 
to the anharmonicity of the trap and therefore oscillates 
slightly more slowly than the parts of the cloud at large 
\z\. Hence, after a few oscillations, the oscillations near 
z = and at large \z\ become out of phase, which results 
in a bending of the cloud in the x-z plane, as illustrated 
in Fig. [5jb). A suitable measure for such a deformation 
is the expectation value of the operator 



<-lb<; 



(Z 2 - (Z 2 ) eq )3 



(33) 



Like qdip and q oct , this operator is defined in such a way 
that it is insensitive to a translation of the cloud as a 
whole in x direction. 

Within the first-order ansatz ([23]) . it is easy to see that 
in a harmonic trap without mean field and to linear order 
in the perturbation, the density oscillates as p(x, y, z, t) — 
p eq {x- (x)(t),y,z), with (x)(t) = ci(t)T/(mujl), i.e., the 
cloud shape remains unchanged. However, in the pres- 
ence of a mean field U, the cloud starts to deform during 
the oscillation. For illustration, the measures discussed 
above for the cloud skewness, (qdip), and for the bending 
(ibend), obtained with the first-order ansatz in the case of 
a harmonic trap (Vo — > oo), are displayed in Fig. [6ja) by 
the dash-dotted and dotted lines. The existence of these 
distortions is in contradiction to the Kohn theorem and 
an unphysical consequence of the crude ansatz. At third 
order, owing to the additional degrees of freedom of the 
ansatz (|2"51) , the unphysical change of the cloud shape in 
the harmonic potential is strongly suppressed, as shown 
by the solid and dashed lines in Fig. [UJa). 

While unphysical in a harmonic trap, a change of the 
cloud shape during the sloshing motion is expected in 
an anharmonic trap. From the preceding discussion it is 
clear that at least the third-order ansatz is required to 
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FIG. 6: Skewness — Im^ip) and bending — \m(qi, en a) of 
the cloud during the sloshing motion excited by V = ax 
at T/Tf = 0.4. In order to have comparable orders of 
magnitude, we multiplied the bending strengths by a factor 
(lj z /uj x ) 2 . (a) In a harmonic trap, the unphysical skewness 
and bending of the first-order approximation (dash-dotted 
and dotted lines) are strongly suppressed when the third or- 
der moments are included (solid and dashed lines). Since 
in a harmonic trap there is no damping, we replaced tj by 
uj + i0.001cu± to generate this plot, (b) Skewness and bend- 
ing during the sloshing motion in the anharmonic trap within 
the third-order approximation. 



get a meaningful description of this effect. In Fig. [BJb), 
third-order results for the expectation values of qdi P and 
qbend during the sloshing motion are displayed. We note 
that, compared to the third-order results for a harmonic 
trap, both responses are now significantly increased, so 
that the peaks correspond to a real physical effect. More- 
over, the observed shapes are typical of what one can 
obtain in a schematic model by coupling two damped 
modes (dipole and bending) to an undamped one (slosh- 
ing) . Finally one can also see that the peaks are centered 
around the sloshing frequency, which is lower than in the 
harmonic case. 



VI. COLLECTIVE MODES RELATED TO 
CUBIC PHASE-SPACE MOMENTS 

In the previous section we have discussed coupling of 
the sloshing mode to other collective modes due to the 
trap anharmonicity. In particular, the ansatz (|28[) con- 
tains the necessary terms to describe the transverse com- 
pressional dipole and bending modes corresponding to 
the operators (|3"l"j) and ([3"3"| . In this section, we will 
study these modes in more detail. However, one should 
keep in mind that, since the excitation operators are al- 
ready third-order ones, the ansatz ([28]) can only give the 
leading-order result for these modes. 

Let us start with the radial dipole oscillation excited 
by the operator (f5T|) . In an ideal Fermi gas in an har- 
monic trap, this operator would excite two modes with 
frequencies uj/lu± = 1 and 3. 

In Fig. the corresponding response function is shown 
for various temperatures. We see that in the collisionless 
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FIG. 7: Strength function — Im{qdi p )(u}) for the radial dipole 
mode (excitation V = aqdip) in the anharmonic trap for vari- 
ous temperatures. The parameters are the same as in Fig. [3] 



regime, i.e., at high temperature, both modes exist but 
their frequencies are significantly lowered by the anhar- 
monicity of the trap. Even at T/Tp = 1.2 both modes are 
still very strongly damped (damping rate ~ 0.3wj_). At 
lower temperature, the damping gets so strong that both 
modes disappear completely: the hydrodynamic regime 
is not reached at any temperature. By artificially increas- 
ing the collision cross section to reach the hydrodynamic 
regime at T/Tp — 0.4, we found that there would be 
a mode at ~ 2uj±_ (~ 2.4wj_ in the case of a harmonic 
trap), but in order to reach this regime in an experi- 
ment one would have to use a trap with a much lower 
radial frequency (i.e., a more spherical or even pancake- 
shaped trap). At the frequency of the sloshing mode, 
the response exhibits a characteristic dip due to the cou- 
pling between the dipole mode and the almost undamped 
sloshing mode (the sloshing is damped mainly because of 
this coupling). 

In order to get a better understanding of the character 
of the two modes in the collisionless regime, we display 
their velocity fields in Fig. [8j Although both modes are 
excited by the same operator, their velocity fields are 
quite different. The velocity field of the low-lying mode 
{to « 0.8w_l at T = T F ), cf. Fig. Ufa), confirms that 
this mode comprises both center-of-mass and skewness 
oscillations. The velocity field of the high-lying mode 
(cj s» 2.3lu± at T — Tp), shown in Fig. [Sfb), resembles 
that of the dipole compression modes in atomic nuclei, 
shown e.g., in Fig. 6 of Ref. f25j. 

Although in an axially symmetric trap the octupole 
mode does not couple to the sloshing mode, let us for the 
sake of completeness briefly discuss this mode. Like the 
dipole operator (|3ip . the octupole operator (|32|) would 
excite two modes at uj /u± = 1 and 3 if the gas was colli- 
sionless and in a harmonic trap. The octupole response 
function for the realistic case is displayed in Fig. [§!&). 
As in the dipole response, both modes are considerably 
shifted downwards because of the trap anharmonicity. 
The damping is even stronger than in the case of the 
dipole modes. Both octupole modes have a similar veloc- 
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FIG. 9: (a) Strength function — Im{q oc t)(Lo) for the radial 
octupole mode (excitation V = aq ct) in the anharmonic trap 
for various temperatures. The parameters are the same as 
in Fig. [3] (b) Velocity field corresponding to the high-lying 
mode at T — Tf (see caption of Fig. [8] for details). 



ity field, as an example Fig.^b) shows that of the higher- 
lying mode. In a nuclear physics context it was found 
that their difference lies mainly in the non-diagonal pres- 
sure tensor (quadrupole moments in momentum space) 
26]. At T/T F = 0.4, there is only a very broad peak 
around the hydrodynamic frequency uj/uj± = \/3 pre- 
dicted in Ref. [27| for a harmonic trap. However the 
strong damping shows that one is still far from the hy- 
drodynamic regime. By increasing again artificially the 
cross-section, the hydrodynamic frequency is found to be 
lowered to ~ 1.58a; j_ due to anharmonicity. 

Finally, in Fig. [TU] we display the strength function 
and the velocity field of the bending mode. In a non- 
interacting gas in a harmonic trap, there would be three 
modes at w = uj x and u> = to x ± 2w z . In the realistic case, 
these peaks cannot be resolved and one sees a damped 
oscillation with a frequency close to u> x . In fact, the be- 
havior of its frequency and damping rate is qualitatively 
similar to that of the sloshing mode, but the temperature 
dependence of the frequency is stronger. Above 0.5Tp, 
the damping is weaker than that of the sloshing mode, 
but this is probably an artefact of our approximation: 
We expect that the bending mode will receive an addi- 
tional damping (analogous to that of the sloshing mode) 
from its coupling to a radial dipole mode modulated in 
z [operator (z 2 — (z 2 ) e q)ldip\i which is not included in 
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the third-order ansatz. Except at the lowest tempera- 
ture, the mixing between sloshing and bending mode is 
weaker than that between sloshing and dipole mode, and 
the bending response does not show the characteristic dip 
that was found in the dipole response. The velocity field, 
shown in Fig. [TOtR ) confirms that this mode corresponds 
to a bending of the elongated cloud in the x-z plane. 

VII. CONCLUSIONS 

In this paper, we studied the frequency shift and 
damping of the sloshing mode of a normal-fluid Fermi 
gas in a realistic anharmonic trap potential. We used 
the moments method in order to find approximate so- 
lutions of the Boltzmann equation with mean field and 
in-medium cross-section. Already at first order, the mo- 
ments method predicts a downwards shift of the fre- 
quency of the sloshing mode. However, the first order 
ansatz is insufficient to describe the coupling between 
the center-of-mass and the internal degrees of freedom 
of the gas, which is responsible for the damping of the 
sloshing mode in an anharmonic trap. We therefore ex- 
tended the ansatz to include also third-order phase-space 
moments. The third-order ansatz contains not only the 
sloshing mode, but also more complicated radial dipole 
and bending modes. The important mechanism for the 
damping of the sloshing mode is its coupling to these 
modes (especially to the dipole mode) due to the trap 
anharmonicity. In addition, the third-order ansatz al- 
lowed us to discuss the radial octupole mode. 

The comparison with the available experimental data 
for the sloshing mode shows that the essential features of 
the frequency shift can be reproduced. Quantitatively, 
the calculated ratios u) s i os h/w x are somewhat lower than 
those given in Ref. Q. Since both the mean field and 
the third-order moments give only small corrections to 
the frequency shift predicted by the first-order moments 
method, further improvements of the theory will prob- 
ably give even smaller corrections. An effect which has 
not been considered in our linear-response study is the 
dependence of the sloshing frequency on the amplitude of 
the mode. However, this effect goes into the wrong direc- 
tion, since the restoring force gets weaker with increasing 
amplitude (see Fig. [T]) and therefore the frequency of the 
sloshing mode will be further reduced (cf. also [1J|)- As 
discussed in Sec. IV Al a possible explanation of the dis- 
crepancy is that the experimental trap frequency ui x may 
have been slightly underestimated in the analysis of Ref. 
0- Actually, our results could be used to determine the 
true value of w x . 

The agreement between the calculated and measured 
damping rates of the sloshing mode is satisfactory and 



makes us confident that our approach contains the essen- 
tial physics to describe the coupling between the slosh- 
ing and internal degrees of freedom. This may also help 
to better understand other effects where the dissipation 
of center-of-mass kinetic energy plays a role, such as the 
heating of the cloud due to laser-beam-pointing noise dis- 
cussed in Ref. jH[. 

In addition to the sloshing mode, we studied the ra- 
dial dipole and octupole modes and the bending mode 
in the x-z plane. However, one should keep in mind 
that for these modes the third-order ansatz contains only 
the leading order, and higher-order corrections can be 
important. We found that, with the present trap pa- 
rameters, the dipole and octupole modes do not be- 
have hydrodynamically at any temperature. The fre- 
quencies of the dipole, octupole, and bending modes are 
more strongly affected by the anharmonicity than the 
frequency of the sloshing mode. For the dipole and oc- 
tupole modes, this can be intuitively understood, since 
these higher-order modes are more sensitive to what hap- 
pens at larger r± , where the atoms feel the anharmonic- 
ity of the trap. The same argument should also apply to 
the radial quadrupole and scissors modes not studied in 
the present paper. One should therefore be careful when 
comparing theoretical results obtained for a harmonic 
trap with experimental ones (as done in 0,111 El), since 
the anharmonicity effects cannot be completely absorbed 
in a renormalized radial trap frequency. 

As a final remark we note that a similar third-order 
ansatz as the one used here [Eq. (|28|) with (x,p x ) <H- 
(z,p z )] can be used to describe one of the higher-order 
axial modes studied recently by the Innsbruck group [29[ 
(the k — 2 mode in the notation of |29(, corresponding 
to a motion as shown in Fig. |4jb) but in z instead of 
x direction). Having more nodes than the usually con- 
sidered collective modes, these kinds of modes are very 
interesting since they are closer to sound waves in uni- 
form systems. 
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